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Abstract 

In this paper we find an exact analytical expression for the number of span- 
ning trees in Apollonian networks. This parameter can be related to signif- 
icant topological and dynamic properties of the networks, including perco- 
lation, epidemic spreading, synchronization, and random walks. As Apol- 
lonian networks constitute an interesting family of maximal planar graphs 
which are simultaneously small-world, scale-free, Euclidean and space filling 
and highly clustered, the study of their spanning trees is of particular rel- 
evance. Our results allow also the calculation of the spanning tree entropy 
of Apollonian networks, which then we compare with those of other graphs 
with the same average degree. 

Keywords: Apollonian networks, small-world graphs, complex networks, 
self-similar, maximally planar, scale-free 



1. Introduction 

In the process known as Apollonian packing [10], which dates back to 
Apollonius of Perga (c262-cl90 BC), we start with three mutually tangent 
circles, and draw their inner Soddy circle (tangent to the three circles). Next 
we draw the inner Soddy circles of this circle with each pair of the original 
three, and the process is iterated, see Fig. [TJ An Apollonian packing can be 
used to design a graph, when each circle is associated to a vertex of the graph 
and vertices are connected if their corresponding circles are tangent. This 
graph, known as Apollonian graph or two-dimensional Apollonian network, 
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was introduced by Andrade et al. [l| and independently proposed by Doye 
and Massen in 0]. 

Apollonian networks have received much attention in recent years as 
two-dimensional Apollonian networks are scale-free, small-world, maximally 
planar, Euclidean and space filling Dynamical processes taking place 

on these networks, such as percolation, epidemic spreading, synchronization 
and random walks, have been also investigated, see [30, Ea, 0, U, 31 1. 
Some authors even suggest that the topological and dynamical properties of 
Apollonian networks are characteristic of neuronal networks as in the brain 



cortex 
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In this paper, we study the number of spanning trees of two-dimensional 
Apollonian networks. This study is relevant given the importance of the 
graphs, and because the number of spanning trees of a finite graph is a 
graph invariant which characterizes the reliability of a network 0j and is 



related to its optimal synchronization and the study of random walks [14]. 
The number of spanning trees of a graph can be obtained from the product of 
all nonzero eigenvalues of the Laplacian matrix of the graph 0] (Kirchhoff 's 
matrix-tree theorem). However, although this result can be applied to any 
graph, this calculation is analytically and computationally demanding. 

Here, we follow a similar approach than for the Sierpihski gasket and 
Hanoi graphs fl, M, EH, 0] and our method to compute the number of 
spanning trees in Apollonian networks is based on the self-similarity of the 
graphs. 



2. Apollonian networks 

We recall in this section the definition and main topological properties of 
two dimensional Apollonian networks. We use standard graph terminology 
and the words "network" and "graph" indistinctly. 

Definition 2.1. An Apollonian network A(n), n> 0, is a graph constructed 
as follows: 

For t = 0, A(0) is the complete graph K% (also called a 3-clique or 
triangle). 

For n > 1, A(n) is obtained from A(n — 1): For each of the existing 
subgraphs of A(n — 1) that is isomorphic to a 3-clique and created at step 
n — 1, a new vertex is introduced and connected to all the vertices of this 
subgraph. Figure shows this construction process. 

We note that Apollonian networks belong to the family of graphs known 
as k-tvees, introduced in the 1970s 0], and which have been extensively 
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n = 2 n = 3 

Figure 2: Apollonian graphs A(n) produced at iterations n = 0, 1, 2, 3 and 4. 
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studied since. They are denned as follows: a /c-clique is a fc-tree; if G is a fc- 
tree, adding a vertex to G and joining it to all vertices of an existing fc-clique 
in G yields a new graph G' , which is also a fc-tree. Thus, Apollonian networks 
are fc-trees, since they start with a 3-clique (that is, the only possible fc-tree 
with 3 vertices), and since each time a vertex is added, it is joined to all the 
vertices of an existing 3-clique. fc-trees are of interest as they are maximal 
graphs for a given tree-width. 

Thanks to the deterministic nature of the graphs A(n), it is possible to 
give exact values for relevant topological properties of these graphs, includ- 
ing degree distribution, degree correlations, clustering, diameter and average 
distance. 



Order and size of A(n).- The order and size of an Apollonian graph A{n 
are V n = \V(n)\ = §(3" + 5) and E n = \E(n)\ = |(3 n + 1). 

Planarity- A graph is planar if it can be drawn on the plane with no 
edges crossing. A planar graph is maximal, or maximally planar, if it cannot 
be extended to a larger planar graph by adding an edge. 

Apollonian networks are maximally planar [3] and this is an important 
feature for its relevance in the development of efficient algorithms [3J. 

Degree distribution. - The cumulative degree distribution of an Apollo- 
nian network A(n) is P(k) oc A; 1-7 , with 7 = In 3/ In 2 = 2.58496, see 
[3,0]. The graph is scale-free. Many real networks share this property with 
exponent values in the same range as A(n) [lj|. 

Correlation coefficient. - From the Pearson correlation coefficient for the 
degrees of the endvertices of the edges of A(n), Doye and Massen ob- 
tained in [3] the exact value of the correlation coefficient and it is always 
negative and thus the network is disassortative. For large values of n, 
r(n) ~ — tt (|) • Thus, this parameter goes to zero as the order of the 
graph increases. Technological and biological networks are usually disassor- 
tative as it is also the case of some information networks, see 
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Modularity and self -similarity.- From Fig. [2J we see that Apollonian 
networks are self-similar, suggesting an alternative way to construct them. 
As shown in Fig. [3j A{n + 1) can be obtained by joining three replicas of 
A(n), labeled by A\, A 2 n and A^, and identifying three pairs of edges. Thus, 
the graphs are modular. Modularity can be quantified with the function 



Q introduced by Newman and Girvan [16]]. For Apollonian networks Q 
0.5938. Moreover, the communities are spatially localized [3, 0]. 
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Figure 3: Recursive construction of Apollonian networks, pointing out their self-similarity. 
A(n + 1) can be obtained by joining three replicas of A(n), labeled here A\, A\ and A\, 
after merging three pairs of edges. 



Clustering.- The clustering coefficient C of A(n) is 0.8284 for large n. 
Furthermore, the clustering coefficient of a vertex of A(n), for large n, is 
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Q| 1 4(fe - §j 
. 1281] : C(k) = —— r-. This property 

is considered to be a signature of a hierarchical network structure, see [2| 

Average distance and diameter.- The modular recursive construction of 



A{n) allows to obtain the exact analytical value for its average distance 24]. 
The average distance of A(n), for n large, follows d(n) ~ In \ V n \ which shows 
a logarithmic scaling with the order of the graph. As the diameter has a 



similar behavior 23(, the graph is small- world. 



3. The number of spanning trees in Apollonian networks 

In this section, we find the number of spanning trees of the Apollonian 
network A(n). For this calculation we use a decimation method [12j which 
has been used to find the number of spanning trees in other recursive graph 
families like the Sierpihski gasket 0], the pseudofractal web 26], and 



some 



fractal lattices 271 ] . The particular structure of Apollonian networks, allow 
us to write recursive equations for the number of spanning trees which are 
solved by induction. 

In the following, we denote by V n and E n the number of vertices and 
edges of A(n). Then a spanning subgraph of A(n) is a subgraph with the 
same vertex set as A(n) and a number of edges E' n such that E' n < E n . 
A spanning tree of A(n) is a spanning subgraph which is a tree and thus 

K = K - 1. 

We call "hub vertices" the three outmost vertices in the construction as 
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shown in Fig. [2] and "hub edges" the three exterior edges which connect the 
hub vertices. 

To simplify our calculations, we define some quantities as follows: We 
denote as s n the total number of spanning trees of A(n). a n represents the 
number of spanning subgraphs of A(n) which consist of three trees such 
that each hub vertex of A(n) belongs to a different tree. In a similar way, b n 
will denote the number of spanning subgraphs of A{n) which consist of two 
trees such that no hub edges belong to the spanning subgraph and one of 
the hub vertices of the subgraph belongs to one tree and the other two hub 
vertices are in the second tree. In a similar way, and taking into account 
the tree to which a given hub vertex belongs, we define b' n and 6". Finally, 
c n is the number of spanning trees of A(n) which have no hub edges. We 
notice that the equation b n = b' n = b" n holds as a result of symmetry, thus, 
in what follows, we will use b n instead of b' n and 6". Figure H] illustrates 
schematically these five types of spanning subgraphs which we denote in all 
figures as A n , B n , B' n , B", C n , respectively. 




Figure 4: Schematic illustration of the five types of spanning subgraphs A n , B n , B' n , B'^ 
and C n derived from A(n). In this figure, vertices represent the hub vertices of A(n) and 
two hub vertices joined by a black line belong to different trees. A red line means that 
the two hub vertices are in the same tree but have no hub edge joining them. 

In the previous definitions, we have not considered the cases where the 
spanning subgraph contains hub edges. We deal with these cases in the 
following three lemmas. 

Lemma 3.1. The number of spanning subgraphs of A{n) which consist of 
two trees such that one hub edge with its two hub vertices belongs to one tree 
while the third hub vertex of A{n) is in the other tree equals a n . 

Proof. Let B n be the set of subgraphs considered in this Lemma. We 
verify the correctness of the result by showing that there exists an one-to- 
one correspondence between the set B n and A n as it is shown in Fig. 
For every spanning subgraph in B n , if we remove the hub edge, then the 
three hub vertices will belong to three different trees, so it belongs to A n . 
Conversely, for every spanning subgraph in A n , if we add a hub edge, then its 
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Figure 5: Graphical illustration of Lemma 13. II The green line means that the hub edge 
joining these two hub vertices belongs to the spanning subgraph. 

two hub vertices belong to one tree and the subgraph is in B n . Thus, there 
exists a one-to-one correspondence between B n and A n , and the cardinality 
of B n is a n . □ 

Lemma 3.2. The number of spanning subgraphs of A(n) such that they 
contain just one hub edge and one hub vertex which is connected to one of 
the hub vertices of this edge through edges of the tree is b n . 

Proof. Let C n be the set of subgraphs considered in this Lemma. As 



Figure 6: Graphical illustration of Lemma 13.21 The green line means that the hub edge 
joining these two hub vertices belongs to the spanning subgraph. 

before, we prove this Lemma by showing that sets C n and have an one- 
to-one correspondence. This can be verified by deleting the hub edge, see 
Fig. [6l Thus, the cardinality of C n is b" n . □ 

Lemma 3.3. The number of spanning subgraphs of A(n) that include two 
hub edges is a n . 

Proof. Consider the bijection between C' n , the set of subgraphs which 
contain two hub edges, and A n (Fig. [7J. □ 

Next we establish a recursive relationship among the parameters s n , a n , 
b n , and c n . 
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Figure 7: Graphical illustration of Lemma 13.31 The green line means that the hub edge 
joining these two hub vertices belongs to the spanning subgraph . 

Lemma 3.4. For n > 0, a n+ \ = 3a„ + 6a^6 n . 

Proof. We prove this result by considering a graphical version of the equa- 
tion (Fig. [8J which represents the recursive construction method of A(n + 1) 
from A(n) and enumerates all possible contributions to a n +i- 

In this representation we only draw four vertices in each case, since each 
non drawn (interior) vertex connects at least to one of these four vertices 
(although they do not have necessarily to be adjacent). This is sufficient to 
determine whether each case belongs to A n+ i, B n+ \ or C n+ \ . 




Figure 8: Configurations needed to find a n +\ from a n and b„. In this representation, the 
red curve denotes the spanning tree to which the two hub vertices belong to. The number 
on the right of the figure counts configurations that, by symmetry, contribute to a n +i 
through the merging process 

Next we should prove that each configuration is correct, but we only 
analyze in detail the first additive term as the other term can be verified 
in a similar way. For this case (see Fig. [9j), hub vertices hi and /14 are 
connected while /12 and /13 are not. Thus, after merging, there are three 
spanning trees and we confirm that this subgraph belongs to A n+ \. Because 
of the symmetry, /14 can also be connected to /12 or /13 and we count three 
times this case. □ 

Lemma 3.5. For n > 0, b n+ \ = + 7a^b n + 7a n b 2 n + a^c n . 
Proof. We prove the lemma by enumeration. Figure [10] shows all the 
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Figure 9: The first configuration which contributes to a„+i. 




Figure 10: Configurations needed to find b n +i- 

distinct possibilities. Again, we only analyze the first case. We label the 
four hub vertices in the same way as in Fig. [9j In the first case, hi, hi, h± 
are all connected while /13 is not. There are two spanning trees, and one 
has no hub edges, so this configuration belongs to set B n+ \. Symmetries 
generate equivalent configurations and the factor is one. □ 

Lemma 3.6. For n > 0, c n+i = + Ylc? n b n + 36a n 6^ + 146^ + 3o^c n + 
12a n 

Proof. As in former lemmas, the proof is by enumeration of all possible 
contributions to c n +i, see in Fig. [TT] the details. In the first case, hi, hi, h§ 
and hi are all connected and the merging process produces a spanning tree. 
As no hub edges are included in it, we can see that this case belongs to set 
C n +i. Besides, because of the symmetry, only this configuration is relevant. 
All other cases are analyzed similarly and we omit the details. □ 

Lemma 3.7. For n > 0, s n+ i = 16a^ + 72a? n b n + 7Sa n b 2 n + 146^ + Q a 2 n c n + 
12a n 

Proof. Figure [121 shows all the configurations contributing s n +i. We do 
not give the calculation details as they are like in former lemmas. □ 
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Figure 11: Configurations that contribute to c n +i. 




Figure 12: Illustration for the calculation of s n +i- Here, a blue line indicates that the 
corresponding hub vertices belong to the same spanning tree (including both the cases 
where they are adjacent and they are connected through other vertices). 
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Lemma 3.8. For n > 0, a n c n — 3b±. 

Proof. By induction. For n = 0, using the initial conditions ao = 1, &o = 
and Co = 0, the result is true. 

Then we assume that for n = k, the lemma is true. For n = k + 1, and 
using Lemmas I3.4H3.6|, we have that 



a-k+iCk+i — 3b k+1 



{3a\ + 6a 2 k b k )(a 3 k + 12a 2 k b k + 36a k b 2 k + Ub\ + 3a|c fc 
+12a fc 6 fc Cfc) - 3(a| + 7a|o fc + 7a fc &| + a 2 k c k ) 2 
3a 2 k (a 2 k + Aa k b k + 7b k - a k c k )(a k c k - 3b\) , 



As by induction hypothesis, we have that a k c k — 3b 2 = we obtain the 
result. □ 



Lemma 3.9. For n > 0, 



O.n+1 



5" 



Proof. Define f n = ""j -1 . From Lemma 13.4^ we have 
«n+i _ 3a^ + Ga^bn 



3 + 6^, 



and thus 

We use Lemmas 
b n +i 



K 
a, 



and 



1 / an+i 
6 V a l 
to obtain: 



1 + 7— + 10 



1 + 7 

an+i 
" 2a3 



1 / a n+ i 
6 



+ 10 



1 f a n +i 
6 



+ 



5a; 



n+l 



18fl6 ' 



which gives: 
b n +l 



a n +i 



bn+l a n 
3 



a n+ i 5a, 



2a3 



+ 



2 

n+l 



18a° 



Thus 



a n +2 



3 + 6 



Vt-1 
a n +i 



3 + 6 



a n +i 
5a n+ i 



1 , 5a n+ i 
"2 18a3 



1 

2 18a2 



5a 



n+l 



3a3 
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which we can write as 

fn+l = 

From the initial condition /o = % = 3 = ■Mr, we finally find 



fn 



3" 

a n+ i 5 n 



□ 



1 , 3' | n 1 , 3" n 

Lemma 3.10. For n > 0, a n , = 3 4+ 4 + 25 4" l_ 4 2. 

Proof. From the former Lemma we have a n +i = 3 n_i a 3 , which with initial 
condition ao = 1 gives 

_1 1 3^_ 1 n 1 I 3" n 
a n = 3 * * ^2 5 4^4 2 . 

□ 

Lemma 3.11. For n>0,b n = ^l^- 1+yi ~ 2n) (5 n - 3 ra ). 

a —3a 3 

Proof. From LemmalJlwe know that b n = ""^ Mn and using Lemma r3.1UI 
we obtain 



6 n = Al5l(-l+3' l -2n) (5 n_3n ) _ 



l 1 

On = 

□ 

Lemma 3.12. For n > 0, c n = i33(3+3"-6n) 5 |(-i+3"-2n) _ 5 n)2_ 

2^2 

Proof. According to Lemma 13.81 we have c n = — 7-, and using the results 
of Lemmas 13.101 and 13.114 we obtain 

Cn = I 3 3( 3 + 3n - 6 ™)53(- 1 + 3n - 2 ™) (3™ - h n f . 

□ 

Theorem 3.13. The number of spanning trees of A{n) is 

Sn = l 3 |(-l+3™- 1 -2(n-l)) 5 i(-3+3"-2(n-l))^3n + 5 n)2_ 
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Proof. The result comes from Lemmas 13.71 13.101 13.111 and 13.121 □ 



After having an explicit expression for the number of spanning trees 
of A(n), we can calculate its spanning tree entropy, which is defined as 



m 



mm 



z = lim — — — . 

The spanning tree entropy is an interesting parameter charactering the 
network structure. 

In 15 

Corollary 3.14. The spanning tree entropy of Apollonian networks is — - — . 

Proof. Define z n = ^r 1 . Prom Theorem 13. 13\ we have 

-8 In 2 + In f + 3 n In 15 - 2n In 135 + 8 ln(3 ra + 5 n ) 
Zn ~ 2(5 + 3 n ) ' 

and finally 

In 15 

z = hm z n = . 

n— yoo 2 

□ 

We can compare this asymptotic value of the entropy of the spanning 
trees for Apollonian networks, z = — 1.3540, with that of other graphs 
with the same average degree. For example, the value for the Sierpin- 
ski graph is 1.5694 [5|] and for the 3-dimensional hyper cubic lattice £3 is 
1.6734 18,@]. Thus, the asymptotic value for Apollonian networks is the 
lowest reported for graphs with average degree 6. This reflects the fact that 
the number of spanning trees in A(n), although growing exponentially, do 
it at a lower rate than graphs with the same average degree. 



4. Conclusion 

In this paper we find the number of spanning trees in Apollonian net- 
works by using a method, based on its self-similar structure, which allows 
us to obtain an exact analytical expression for any number of discs. The 
method could be used to further study in this graph, and other self-similar 
graphs, their spanning forests, connected spanning subgraphs, random walks 
and vertex or edges coverings. Knowing the number of spanning trees for 
Apollonian networks allows us to show that their spanning tree entropy is 
lower than in other graphs with the same average degree. 
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